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The topological phases in one-dimensional quantum walk can be classihed by the coin parameters. 

By solving for the general exact solutions of bound states in one-dimensional quantum walk with 
boundaries specihed by different coin parameters, we show that these bound states are Majorana 
modes with quasi-energy E — 0, tt. These modes are qualitatively different for different boundary 
conditions used. For two-boundary system with symmetric boundary conditions, the interaction en¬ 
ergy between two Majorana bound states can be computed, as in the case of a hnite wire. Suggestion 
of observing these modes are provided. 

PACS numbers: 03.65.Ge, 03.65.Vf, 73.63.Nm, 37.10.Jk 


- Introduction Quantum walk provides an excellent 
platform for us to further our understanding of the in¬ 
terference and superposition features that lead to a non- 
classical dynamic evolution. The recent advances in dis¬ 
crete time quantum walk m have created not only 
many interesting research endeavors in its application, 
such as in quantum computation HE], but also in the 
possibility of finding new topological phases in condensed 
matter physics [12]. In discrete quantum walk HHIISl- 
the many parameters characterizing the quantum 
coin allow the theorists to explore different physical sce¬ 
nario for the walker, thereby relating theoretical predic¬ 
tion with experiments for the studies of different dynamic 
scenarios EniiiT]. 

For system with two internal degrees of freedom, the 
quantum coin is a 2 x 2 unitary matrix, where Hadamard 
coin is usually employed. Here we consider a general¬ 
ization of the Hadamard coin with a rotational matrix 
parametrized by a parameter 0 G (—7r,7r), 

cosO sin^ . . 

— sin 0 cos 0 ^ ^ 

The basis for this matrix are which are the 

two internal degree of freedom of the walker, with the 
left (L) and right (R) labels for one-dimensional quantum 
walk. The quantum walk is implemented by an iteration 
of a unitary transformation acting on a two-state particle, 
with a wave function l^ip) = (<^n|^) P bn\R)) \n) defined 

on a discrete one-dimensional space where \n) is the state 
that quantum walk localized at position of integer n, with 
On and bn being the left and right components of the 
wave function at position n. The unitary transformation 
U = SC consists of a quantum coin operation C that acts 
on the internal state, followed by a shift operator S to 
evolve the particle coherently in position space. The shift 
operator S = \L){L\ ^\n){n P 1\ P \R){R\ Z] l^)(^ “ 1| 
moves the left and right component to the left and 
right direction by one unit respectively. The eigen¬ 
states of a quantum walk which coin is translational 
invariant are plane waves characterized by their quasi¬ 


momentum /c, with quasi-energy E and is written as 
'0/c(n, t) = [a/c, bk]^. The quasi-momentum and 

quasi-energy have a periodicity of 27r and we restrict the 
value to be from (— tt, tt) (from now on, we omit the prefix 
’’quasi-” in front of energy and momentum). From the 
evolution equation, the energy E and the momentum k 
satisfies the dispersion relation cosF^ = cos 0 cos k and 
the normalized vector [ok^bk]^ is. 


1 

Y^2 sin E (sin + cos 0 sink) 


i sinOe^^ 

sin E + cos 0 sin k 


( 2 ) 


The dispersion relation shows that for a given momen¬ 
tum /c, there exist two eigenstates with energy E = 
± arccos(cos 0 cos /c), which we denote them by |/c, +) and 
|/c, —) , corresponding to the positive and negative energy 
eigenstates respectively. All the eigenstates form two en¬ 
ergy bands with a gap which closes at = 0 and E = tt 
when 0 = 0 OT 0 = TT. 

- Topological Phases Recently T. Kitagawa has shown 
that quantum walks have a nontrivial winding number of 
1 due to its hidden chiral symmetry and for the split-step 
quantum walks there can be two topological phases [22| . 
Firstly, we write the unitary transformation U{k) for a 
specific momentum k as 


U{k) = 


cosOe^^ sinOe^^ 
— sm0e~'^^ cos0e~'^^ 


( 3 ) 


which can be rewritten in the form of the evolution opera¬ 
tor U{k) = provided that H{k) = E(k) [n(k) • a] 

with E{k) chosen to be positive and the unit vector n(k) 
defined by. 


n(k) = —(sin 0 sin /c, sin 0 cos /c, cos 0 sin k )/ sin E (4) 

One can verify that U{k) = {cosE)I — isinE [n(k) • a]. 
Here I is the identity matrix and a the Pauli matrix. 
The Hamiltonian H{k) can also be written in terms of 
the eigenstates, H{k) = E\k,P){k,P\ — E\k,—){k,—\. 
This Hamiltonian H{k) can be mapped to a unit sphere 
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FIG. 1. The trajectory of vectors n{k) (black circle) and 
A(^) on the unit sphere (blue semi-circle). The vector A(^) 
is perpendicular to the black curves and dehnes the winding 
direction of the vector n{k) 


through the vector n{k) that traces out a circle on the 
equator of the unit sphere when k goes from —tt to tt 
( see Fig. 1) and the circle provides a winding number 
corresponding to the topological number of the quantum 
walk. The sense of the winding of the circle around the 
equator can be defined through the normal vector A{k) 
such that A • n{k) = 0, 

( 5 ) 

As k changes, n{k) moves on a plane with normal A{k). 
This is due to the hidden chiral symmetry in the Hamil¬ 
tonian: 




where H is the chiral symmetry operator which rotates 
n{k) by tt around the axis defined by A{k). This chi¬ 
ral symmetry is essential for the existence of topolog¬ 
ical phase in one-dimension and the resultant topolog¬ 
ical phases are called symmetry protected topological 
phases. With this chiral symmetry, the Hamiltonian can 
be off-diagonalized and the topological number is de¬ 
termined by the off-diagonal terms m, which here is 

achieved by the application of a unitary transformation 
A = g-icr^(7r/2-0)/2^ 


A-^HA = 


E 

sin£^ 


0 h{k) 

h'^{k) 0 


h = sin k — i(sin 0 cos k) 


( 7 ) 

The unitary transformation A is a rotation around y-axis 
(which is along n(0)) by 7r/2 —It transforms H{k) to a 
new Hamiltonian that has only (7^ and Gy components, so 
that as k changes, it traces out a circle on a 2-dimensional 
x — y plane. The topological number can then be defined 
as the winding number of h around the origin. 


The discontinuity of topological number m at 6> = 0,7r 
can be understood as the discontinuity of the vector 
A{0). As 0 increases from 0 to tt, A{0) goes from point 
I to point H through the blue dash semicircle shown in 
Fig. 1. In this regime, A{0) = (cos6>, 0, — sin6>). As 0 
increases from tt to 27r, (or equivalently from — tt to 0) 
A{0) = (—cos6>, 0, sin6>) and we see that it jumps from 
point H to point I and goes through the blue dash semi¬ 
circle once again. To handle this discontinuity, we intro¬ 
duce another continuous vector B(^) = (cos^, 0, — sin^), 
which does not contain the sign function sgn(sin^). For 
0 between 0 and tt, B(6>) = A{0); while for 0 between 
TT and 27r, B(6>) = —A{0). The positive winding orienta¬ 
tion is defined to be along vector B(6>) = (cos 6>, 0, — sin 6>) 
which goes through a complete circle continuously from 
point I to H and then from H to I again. Accordingly, the 
topological numbers should he m = A • B = sgn(sin^). 
The ambiguity in the definition of topological number 
occurs only at 0 = 0,7r because the band gap closes at 
these points. The absolute value of the topological num¬ 
ber is the number of times that the loop winds around 
the equator, which gives a winding number \m\ = 1, but 
there exist two topological phases with topological num¬ 
ber m = ±1. 

- Majorana Boundary Mode From bulk-boundary cor¬ 
respondence [24] , there should be two gapless states local¬ 
ized at the surface when these two topologically distinct 
phases are joined together. These gapless states occurs 
only at E = 0, tt, thus there should be two bound states 
of energy either dit E = 0 or E = n on the boundary 
of quantum walk with sin^ of different sign. Based on 
these observations, we seek to find these boundary states 
by considering a boundary at the origin (n = 0) where 
the quantum walk has different coins on the two sides 
of the boundary. Specifically, 0i is used for n < 0 and 
O 2 0i for n > 1. Following the analysis of Ramsauer 
effect in quantum walk [25]. we look for bound states for 
this quantum walk by the following ansatz, with the con¬ 
dition that the wave function decays to zero at infinity. 


tp{n) 


re 


Kin 


^—iKi 

^—iKi 




^iK,2 

^iK2 


n <0 

n> 1 


( 9 ) 


Here we replace the real k in [a/e, bk]^ with an imaginary 
one. For n < 0, the spinor is determined by Oi] for n > 
1, the spinor is determined by O 2 . Using the dispersion 
relation cosE = cos6>cos/c, we can relate the energy of 
this hypothetical bound state with the decay constant ni 
and 1 ^ 2 ^ 


cos E = cos 0i cosh = cos O 2 cosh 1^2 (10) 


[ dk-^lnh(k) = sgnisinO) 
27ri J_^ dk ^ ^ ^ ^ ^ 


By imposing the continuity conditions on the bound state 
wave function at n = 0,1, we get the relations between 
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the r and t: 


re — te (H) 


where the first equation conies from left component of 
' 0 ( 0 ) and the second comes from the right component of 
0(1). In order for the bound state to exist, we need the 
existence of nontrivial solution of r and t, which needs 
sin 0i and sin O 2 to be of opposite sign (refers to Appendix 
I in the Supplementary Material). This implies that two 
sides of the boundary are topologically distinct. Under 
this condition, there are always two bound states solution 
with energy ^ = 0, tt, which comes in pair. This is also 
predicted from the bulk-boundary correspondence. 

For a concrete example of quantum walk where the 
coin parameter is positive on the left and negative on the 
right of the origin, the bound states solutions are. 


E = 0 0 (n) = 


Xi 

-1 

X2 

-1 


E = IT 'ip{n) 


(-xi)"-i 

(-X 2 )"-' 


Xi 

1 

X2 

1 


n < 0 

n > 1 

n < 0 

n > 1 


( 12 ) 


where xi is ( 1 +sin^i)/cos^i and X 2 is (l+sin^ 2 )/ cos^ 2 - 
The two sides of the wave function only depend on the 
coin parameter of their own sides. For two regions with 
same topological number, the wave function on both re¬ 
gions can be joined together and the continuity condition 
on the boundary will be satisfied automatically. This ob¬ 
servation is important in our discussion for finite wire. 

Now that our ansatz for the boundary state is possible, 
we claim that it is a Majorana bound state, which is a 
0 + 1 dimensional particle whose antiparticle is itself [26] , 
(0 + 1 stands for 0 dimension in space and 1 dimension 
in time). The system has an intrinsic particle-hole sym¬ 
metry. We denote the complex-conjugation operator as 
K. Then, the eigenstate with momentum k and energy 
E and the eigenstate with momentum —k and energy 
—E are related by an antiunitary operator F = —K, 
\k,E) = r| — k,—E). This antiunitary operator is the 
particle-hole operator so that our effective Hamiltonian 
has an intrinsic particle-hole symmetry, 

rH{k)r-^ = {-K)H{k){-K) = -H{-k) (13) 

This particle-hole symmetry shows that a creation of a 
particle with energy E is equivalent to the annihilation of 
a hole with energy —E. Because of the periodicity of the 
energy, E and —E coincide when E = 0, tt, the creation 
and annihilation operator are the same for = 0 and 
E = 7 T and the boundary modes are Majorana modes. 

- Finite Wire We now consider quantum walk on a fi¬ 
nite wire. Note that when 0i = ± 7 r/ 2 , the decay constant 


Ki in the bulk is infinity according to the dispersion re¬ 
lation in Eq. and thus the quantum walk does not 
penetrate through a site with coin parameter 0i = ±7r/2 
which corresponds to 


C = ± 


0 1 
-1 0 


(14) 


This coin interchanges the left and right components with 
or without an extra sign so that it reflects the quantum 
walk. Thus, the proper boundary conditions describing 
the finite wire is to set the coin parameters at the two 
ends of the wire to be ±7r/2, and this resembles the quan¬ 
tum wire model proposed by Kitaev [27|. 

Now consider the general case of a two-boundary sys¬ 
tem: the coin parameter is 0 i for n < 0 , and O 2 for sites 
from n = 0 to n = A/", and ^3 for n > N. We can now 
have two topologically distinct configurations. The sym¬ 
metric configuration is defined by 63 = 0 i as shown in 
Fig.2 (a) and the antisymmetric configuration is defined 
by ^3 = as shown in Fig.2(b). We can reduce this 
general case to the finite wire by letting 0 i approaches 
± 7 r/ 2 . Let’s now use the following ansatz for the bound 
state wave function. 


0 (n) = Ae-^^^ 


= Ce-'^^^ 


^iK,2 




^ — iK2 
^ —m 2 


Q < n < N 
n > N 


= De'^^^ 


^—iKi 


n <0 


(15) 

where is positive so that the wave function vanishes 
at the infinity. The plus sign of the ± sign in the ansatz 
for n > N is responsible for the symmetric configuration 
(Fig.2a), while minus sign is responsible for the antisym¬ 
metric one (Fig.2b). Here, the decay constant and the 
energy are related by Eq. The boundary conditions 
at n = 0 and n = N are 


E Bb—ii^^ — Eb-ii^-^ 

Ae'^^ai^^ + = De~'^^a_i^^ 

(16) 

This system of linear equations has nontrivial solution 
only when its determinant is zero. Thus, we obtain a 
relation between energy and decay constants. Eor the 
symmetric configuration (Eig.2a), this relation is, 

sinh [n 2 {N + l)](sin^ E — sin 6 >i sin 6 ^ 2 ) 

= cosh [/^ 2 (A^ + 1 )] cos 0 i cos O 2 sinh sinh tv 2 


while for the antisymmetric configuration (Eig.2b), it is, 

sinE(cos 6 >i sinh/^i tanh/^ 2 (A^ + 1 ) + cos 6^2 sinh/i: 2 ) = 0 

(18) 
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(a) 


(b) 

FIG. 2. The spatial profile of the coin parameter (top) and 
topological number (bottom) for the symmetric (a) and anti¬ 
symmetric (b) conhguration. The topological prohles for the 
symmetric and antisymmetric conhguration have one and two 
jumps respectively. 


Notice that Eq. is satished when sinE^ = 0, i.e. 


E = 0, TT, however, Eq. 17 does not have solution when 
E = 0, TT. This difference can be understood by looking 
at the spatial prohle of the topological number shown in 
Eig.2. Eor the symmetric conhguration (Eig.2a), there 
are two jumps in the topological number prohle and is 
of special interest. Our discussion on boundary states 
indicates the existence of bound states localized at the 
boundaries with discontinuity in topological number. Eor 
the symmetric conhguration with hnite the system 
has bound states only when sin 0 i and sin O 2 are of oppo¬ 
site sign so that there exist a discontinuity in topological 
number prohle. Eor two-boundary system with symmet¬ 
ric boundary condition, such as the two ends of a wire, 
the two Majorana states at the two ends interact with 
each other and this interaction between them gives rise 
to a small energy difference away from E = 0, tt. We 


substitute Eq. into Eq. and obtain the energy de¬ 
viations from E = 0 , TT of the bound states with different 
separation N and coin parameters O 2 different from 0i = 
1 x12 (see Appendix II of the Supplementary material). 
The small energy deviations decrease exponentially by 
ex.p{—K 2 N) with ^2 = — log |(1 — sgn(^ 2 ) sin^ 2 )/cos^ 2 |- 
This deviation is the interaction energy of the two Ma¬ 
jorana state located at the two boundaries. Eor the anti¬ 
symmetric conhguration (Eig.2b), there is only one jump 
in the topological number prohle and thus, it is equiva¬ 
lent to the topological number prohle of Majorana bound 
state obtained with one boundary. In this case, there are 
always bound state solutions with energy E = 0, tt. Eor 
example, for energy E = 0 when 0i < 0 and O 2 > 0, the 
bound state in the antisymmetric conhguration is. 


A = sin 0i 


B = 0 


C = D = sin 6»2 


(19) 


This bound state solution has B = 0 which means that 



FIG. 3. The probability distribution of boundary modes for 
three diherent conhgurations of a hnite system with length 
A = 10 and O 2 = — 7r/4. The hrst and the second value in 
the legend is the coin parameter on the left and right end 
respectively. 


the bound state is located at n = 0 which is also the 
location of the jump in the topological number prohle 
as shown in Fig. 2b. One can verify the solution with 
E = TT is also located at the discontinuous point of the 
topological number prohle. In Fig. 3 , we illustrate the 
probability distribution of the bound states on a hnite 
wire in three different conhgurations. The coin parame¬ 
ter on the line is O 2 = —7r/4. For the symmetric conhg¬ 
uration which coin parameters are 7r/2 on the two ends, 
the state is localized at the two ends. However, for the 
antisymmetric conhguration which coin parameters are 
7r/2 and —7r/2 at the two ends, the state is localized at 
the end with coin parameter 7r/2. 

- Conclusion We have demonstrated the existence of 
two topological phases in a quantum walk system. We 
then investigated the boundary states for systems with 
single boundary in topological phases. Through an 
ansatz for bound state, we hnd solution for the bound¬ 
ary state and show that they are Majorana states with 
energies E = 0, tt. For systems with two boundaries in 
topological phases, such a wire, we hnd same solution 
as the single boundary system if the boundary condition 
is antisymmetric. However, if the boundary condition is 
symmetric for the two boundary system, we hnd a small 
change in energy away from E = 0, tt, which is a signa¬ 
ture of the interaction energy of the Majorana states. For 
the special case of a hnite wire, the energies of the Majo¬ 
rana states have also been calculated as a function of the 
length of the wire. Possible observation of these bound 
states can be made following the work of Kitagawa et 
al.[28] for photonic quantum walk. Our work may shed 
some insights into the recent experimental work on the 
interaction effects on a Majorana zero mode leaking into 
a quantum dot pa EO]. Experiments along these sys¬ 
tems may lead to the observation of Majorana modes or 
boundary modes localized at the end of a hnite quantum 
walk system either using photon or electrons, thereby ver¬ 
ifying our theoretical calculation. On theoretical front. 
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we like to seek analogy in a generalization of topological 
phases in higher-dimensional quantum walk system, as 
well as using these results for quantum computation. 
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Supplementary Material 

Appendix I. Condition for the existence of bound 
state in single-boundary system 


For the existence of the bound state in single boundary 
system, there must be nontrivial solution of r and t of Eq. 


11 Thus we set the determinant of these linear equation 


to zero and this yields a condition between the energy E 
and the decay constants: 


i sinF^ 


sin $2 cos 0 i sinh ki + sin 0 i cos O 2 sinh k ,2 
sin 0 i — sin O 2 


(A.20) 

As the right hand side of the equation is real while the left 
hand side is imaginary, equality only hold when sinF^ = 
0, implying E = 0 or tt. The decay constant k in this 
case is 


f; = 0 , e'" 

E = 7v, 


1 ± sin 6> 

cos 0 
1 ± sin 0 

cosO 


(A.21) 


Since at infinity the wave function decays to zero, k,i and 
tv 2 both must be positive, so that exp(/i:) pick up the 
plus sign of the ± sign for positive sinO and the minus 
for negative sin 0. Using these /i: for F^ = 0, tt in Eq. |A.2Q[ 
we obtain an equality. 


0 = sin O 2 sin 0 i 


sign(sin^i) + sign(sin^ 2 ) 
sin 0 i — sin O 2 


(A.22) 


However, this equation holds only when sin 0i and sin O 2 
are of different sign. Thus, our hypothetical bound state 
exists whenever the topological phases on the two sides 
of the boundary are different. Furthermore, there are 
always two bound states solution with energy E = 0, tt, 
which comes in pair in agreement with the prediction 
from the bulk-boundary correspondence. 


Appendix II Condition for the existence of bound 
state in two-boundary system in symmetric 
configuration 


For two-boundary system in the symmetric configu¬ 
ration in the limit when N tends to infinity, we have 
sinh/i: 2 Y sign(/^ 2 ) cosh/i: 2 Y and Eq. 17 becomes. 


sin^ E — sin 0i sin O 2 
cos 0 i cos $2 


= sgn(/^ 2 ) sinh sinh tv 2 (A.23) 
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TABLE L The bound state energies E{02) of finite systems 
with length iV, whose coin parameters are O 2 in the middle 
and 0i — —7^12 at the two ends as shown in Fig.2a. The 
system has four bound states solution with energy E, — E, 
TV — E and —(tt — E). 


N 

E(7r/3)/7r 

E(7r/4)/7r 

E(7r/6)/7T 

1 

2.13 X 10“2 

4.68 X 10“2 

8.04 X 10“2 

2 

5.69 X 10“^ 

1.89 X 10“^ 

4.31 X 10“^ 

3 

1.52 X 10“® 

7.77 X 10“® 

2.41 X 10“^ 

4 

4.08 X 10“^ 

3.21 X lO”'^ 

1.37 X 10“2 

5 

1.09 X 10“^ 

1.33 X 10“® 

7.89 X 10-® 

6 

2.93 X 10“® 

5.52 X 10“"^ 

4.54 X 10“® 

7 

7.85 X 10“® 

2.29 X 10“^ 

2.62 X lO"'^ 

8 

2.10 X 10“® 

9.47 X 10“® 

1.51 X lO"'^ 

9 

5.64 X 10“’’ 

3.92 X 10“® 

8.73 X 10“^ 

10 

1.51 X 10“'^ 

1.62 X 10“® 

5.04 X 10“^ 


Equality can hold only when sin 0i and sin O 2 has opposite 
sign because ki is required to be positive from the bound¬ 
ary condition at infinity. This agrees with the topological 
condition in single boundary system. Now, when sin^i 
and sin 6>2 has opposite sign, we find two solutions with 


energy E = O^tt with, 

sinh/^i = sign( 6 >i) tan 6 >i 
sinh tv 2 = tan O 2 


(A.24) 


Similarly, for finite N^ bound state solution can exist 
only when sin^i and sin ^2 are of opposite sign. From 
the dispersion relation Eq. 10, we see that |sinE^| is 


less than both | sin 0 i \ and | sin O 2 \ for a bound state 
solution and thus (sin^ — sin sin ^ 2 ) is negative if 
sin sin ^2 > 0. Note that the finite condition at in¬ 
finity requires sinh/^i > 0 so that equality of Eq. [T7| 
requires (sin^ E — sin 6 >i sin 6 ^ 2 ) > 0. Thus, for finite wire 
bound state solutions can exist only if sin 6 >i sin 6^2 < 0 . 
Using Eq. Eq. can be written as an equation 
independent of /^i. 


sinh [k, 2 {N + l)](sin^ E — sin 6 >i sin 6 ^ 2 ) 

= cosh [k, 2 {N + 1)] cos O 2 sinh ti 2 \/ cos^ E — cos^ 0i 

(A.25) 

This new equation can be also applied to finite system 
where 0i = 7 r /2 and ki = + 00 . The equation can now 
be solved by using Eq. [^to relate ki and E. 







